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Measurements of direct photon production in p+Pb and p+C collisions at ^/snn = 17.4 GeV are 
presented. Upper limits on the direct photon yield as a function of j?t are derived and compared to 
the results for Pb+Pb collisions at ^/snn = 17.3 GeV. The production of the rj meson, which is an 
important input to the direct photon signal extraction, has been determined in the rj — > 27 channel 
for p+C collisions at ^/snn = 17.4 GeV. 
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The search for a deconfined state of nuclear matter, 
the so-called Quark-Gluon Plasma (QGP), has been the 
driving force behind the program of high-energy heavy- 
ion research in the last decades. A variety of results from 
measurements at the CERN-SPS and RHIC at BNL, and 
most recently at the LHC at CERN strongly suggest the 
existence of a deconhned phase. One of the earliest pro- 



posed QGP signatures, the production of thermal direct 
photons [lj, nas been exceedingly difficult to establish 
experimentally 0, 0] ■ 

Since photons do not interact via the strong force, they 
leave the interaction zone without modification, carry- 
ing information about the system at the time of their 
production. Direct photons, i.e. all photons not origi- 



2 



nating from decays of long-lived particles, thus allow to 
probe the entire history of the collision process. The di- 
rect photon yield includes a prompt contributions from 
initial hard scatterings of the partons, that can be de- 
scribed by perturbative QCD calculations. Thermal di- 
rect photons are those produced in the later equilibrated 
phase of the heavy-ion collision with a yield that will be 
characterized by the temperature of the hot and dense 
matter created |4J. Measurements of direct photon pro- 
duction in p+A or p+p collisions are a necessary baseline 
to determine the prompt photon contribution in A+A 
collisions. Aside from nuclear effects on parton distribu- 
tion functions, the prompt direct photon contribution in 
a nucleus-nucleus collision is expected to scale with the 
number of binary nucleon-nucleon collisions, as has been 
demonstrated in recent measurements by the PHENIX 
experiment at RHIC [!, @ . An accurate determination 
of the prompt direct photon component based on a di- 
rect photon measurement in p+p collisions at the same 
incident energy Q, which is in agreement with pQCD 
predictions, allowed to extract a possible excess contribu- 
tion of direct photons in A+A collisions, that may then 
be attributed to a thermally equilibrated phase. The 
spectra of these thermal direct photons can be used to 
constrain the initial temperature of the medium created 
in the heavy-ion collision. Recent measurements of the 
virtual direct photon spectrum in the electron-positron 
channel by the PHENIX experiment have provided first 
evidence for a direct photon enhancement beyond the ex- 
pected prompt photon yield in central Au+Au collisions 
at snn = 200 GeV [7j, which may be explained by a 
thermal production mechanism. 

The first observation of a significant direct photon sig- 
nal in heavy-ion collisions was reported by the WA98 
experiment [1, Q at the SPS. The direct photon yield 
was measured in central Pb+Pb collisions at ^/snn = 
17.3 GeV. At the time, no direct photon measurement 
in p+p collisions at yfs = 17.3 GeV was available as a 
reference to determine the expected prompt photon sig- 
nal. Instead, direct photon measurements in p+p and 
p+A collisions at the higher ^/snn of 19 GeV by other 
experiments were used as a reference to deter- 

mine the prompt photon signal by applying xx-scaling 
(it = 2pr/\/s) to account for the difference in the 
center-of-mass energies [H[l3j]. However, these measure- 
ments were only available for comparisons above pt ~ 2 
GeV/c, and differed from each other by as much as a 
factor of 2. Furthermore, pQCD calculations of prompt 
direct photon production in the SPS energy regime were 
generally unable to reproduce the measurements without 
the introduction of intrinsic fcx contributions 14]. As a 
result, the prompt photon yield could not be estimated 
without large systematic uncertainties, and therefore re- 
liable limits on a thermal contribution to the observed 
direct photon excess in central Pb+Pb collisions could 
not be inferred without similarly large systematic uncer- 
tainty [H[il. Therefore, the magnitude, and even the 
existence, of a possible thermal contribution to the direct 



photon spectrum at SPS energies has remained an open 
question. 

In this paper, results of a direct photon analysis of 
p+Pb and p+C collisions at v / snn = 17.4 GeV by the 
WA98 experiment are presented. Upper limits on the 
direct photon yield are obtained and used in a compar- 
ison to the previously published WA98 results on di- 
rect photon production in central Pb+Pb collisions at 
a/snn = 17-3 GeV. New measurements of the 77-yield 
in p+C collisions are presented and combined with the 
WA98 measurements of the neutral pion yield to extract 
the 77/7r°-ratio, which is a crucial input to the determi- 
nation of the direct photon yield. 

In 1996 the WA98 experiment took data with a sec- 
ondary beam of 160 GeV/c momentum on Pb and C 
targets, corresponding to a nucleon-nucleon center-of- 
mass energy of ^/snn = 17.4 GeV. The primary proton 
beam of 450 GeV/c from the SPS accelerator impinged 
upon a beryllium target and a secondary mixed beam 
of 160 GeV/c momentum, consisting primarily of pro- 
tons and charged pions, was selected and delivered to 
the WA98 experiment. A clean beam trigger was defined 
as a signal in a plastic scintillator start counter, located 
3.3 m upstream of the target, with no coincident signal 
in a veto scintillator counter (which had a 3 mm diame- 
ter circular hole and was located 2.7 m upstream of the 
target), or in beam halo scintillator counters (covering 
the transverse region beyond the veto counter up to a 
radius of 25 cm). Two threshold Cerenkov counters were 
used for the identification of the beam particles. Proton 
projectiles were selected by requiring signals below the 
corresponding thresholds for pions and kaons in both de- 
tectors. Only events where both Cerenkov detectors iden- 
tified the beam particle as below threshold, and therefore 
a proton, were used in this analysis. The minimum bias 
trigger was determined from the total transverse energy 
deposited in the ran ge 3 .5 < 77 < 5.5 as measured by the 
MIRAC calorimeter Q. 

Two targets were used for the data presented here: A 
208 Pb target with a thickness dpb of 0.436 mm and an 
areal density ppt of 495 mg/cm 2 , and a 12 C target with 
d c = 10.022 mm and p c = 1879 mg/cm 2 . A total of 
1.2 • 10 6 (1.0 • 10 6 ) minimum bias events were recorded 
for p+C (p+Pb) collisions, with a factor of ~ 30 (~ 8) 
more minimum bias events sampled by a high energy 
photon trigger described below. The measured WA98 
minimum bias cross section a m b of 193 mb (1422 mb) 
for p+C (p+Pb) corresponds to 86% (81%) of the total 
inelastic cross section. 

Photons were measured with the 10,080 module lead- 
glass electromagnetic calorimeter (LEDA) located at a 
distance of 21.5 m downstream of the target with roughly 
one quarter of full azimuthal coverage in the pseudo- 
rapidity range 2.3 < 77 < 3.0. A Charged Particle Veto 
(CPV) detector was placed 1 m before LEDA to reject 
charged clusters. To enrich highly energetic photons in 
the data sample, a high energy photon (HEP) trigger 
derived from overlapping 4x4 groups of adjacent mod- 
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FIG. 1: Invariant inclusive photon yields from p+Pb and p+C 
collisions at ^/snn = 17.4 GeV. The error bars represent the 
quadratic sum of the statistical and systematic uncertainties. 



ules (each group offset by 2 towers) in the LEDA was 
used. The HEP trigger reached its full efficiency for 
p T > 0.8 GeV/c. With the HEP trigger, an additional 
1.5 • 10 6 (0.5 • 10 6 ) p+C (p+Pb) events were recorded, 
equivalent to an additional 3.9 ■ 10 7 (8.2 ■ 10 6 ) sampled 
minimum bias events. 

This analysis follows the procedure described in 
Ref. The fully corrected inclusive photon spectra 

measured with the electromagnetic calorimeter for both 
the p+Pb and p+C data sets are presented in Fig.Q] For 
Pt > 1-2 GeV/c the triggered data sample is used. 

The spectra have been corrected for the geometrical ac- 
ceptance and efficiency, analogous to previous measure- 
ments by WA98 as described in Refs. [1, [Hj]. The ac- 
ceptance correction accounts for the limited spatial cov- 
erage of the detector for single photons. The acceptance 
is independent of pt and amounts to 0.245 in the rapid- 
ity interval 2.0 < y < 3.2. The efficiency correction has 
been determined from GEANT simulations in which sim- 
ulated clusters have been embedded into the raw data to 
undergo the full reconstruction and analysis procedures 
that are applied to the measured data. The results have 
been validated with an independent Monte-Carlo sim- 
ulation. The efficiency correction accounts for defective 
detector modules, finite energy resolution effects, and the 
cuts applied to the data. The absolute energy scale was 
fixed to the neutral pion mass as described in Ref. [l8j |. 

Off-target events were subtracted on a statistical basis 
by using the results obtained from data samples taken 
without a target in place. The CPV was used to iden- 
tify photon conversions and determine the fraction of 
charged particles registered as photons in the LEDA. 
As the CPV was not available during all p+A runs, the 
fraction of charged particles has been determined on a 
statistical basis and with a limited pt coverage. It has 
been extended by scaling the correction used in periph- 
eral Pb+Pb collisions, the scaling factors have been val- 
idated with Hijing and AMPT simulations. The correc- 
tion for the amount of (anti-) neutrons falsely detected 



in the LEDA as photon hits was determined from sim- 
ulations based on the results for peripheral Pb+Pb col- 
lisions presented in Ref. [||, using scaling factors deter- 
mined from the same Hijing and AMPT simulations that 
were used for the CPV information. To account for the 
steeply falling shape of the spectra, the contents of the 
Pt bins are shifted along the y-axis employing a fit with 
a Hagedorn function fuag = a ■ [b/(b + Pt)] 2 to represent 
the correct value at the bin center. 

The systematic errors on the inclusive photon measure- 
ment are summarized in Table U for two representative 
Pt values. The uncertainty of the efficiency correction 
includes the uncertainties introduced by the cut on the 
shower shape in the detector, by the finite energy resolu- 
tion of the detector, which is the dominant uncertainty at 
the highest pt, as well as the conversion correction. The 
Pt dependent uncertainty for the energy scale has been 
determined by variation of the scale within its ±1.5% ac- 
curacy, estimated from the neutral pion analysis. The 
errors bars of Fig. [1] show the quadratically summed sta- 
tistical and systematic uncertainties 
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FIG. 2: Invariant r\ yield measured in p+C collisions at 
\/ snn = 17.4 GeV in the r\ — > 27 channel. The error bars 
represent the quadratic sum of the statistical and systematic 
uncertainties. 

The direct photon yield is calculated as a fraction of 
the inclusive photon yield presented above, following the 
method described in Ref. Q • The direct photon yield for 
each pt interval is calculated as: 



Tdircct — Tinclusivc Tdecay 
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In the double-ratio i? 7 , systematic uncertainties that 
are present in the neutral pion and in the inclusive pho- 
ton measurement partially cancel, thereby reducing the 
total systematic uncertainty on the measurement. The 
neutral pion yield in the 7r° — > 27 channel for the two 
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data sets has been presented in Ref. [18j ■ We use a fit to 
the 7T° spectra with a Hagedorn function as reference in 
the measured (7/7r°) meas ratio to overcome the different 
bin widths. 

To determine the (7 decay ratio for the decay pho- 
tons a Monte Carlo simulation was employed to calculate 
the photonic decays of the relevant mesons - 7r°, ?y, ?/, w, 
and K® - based on the measured neutral pion spectra in 
p+Pb (p+C) collisions and by employing mx-scaling for 
the heavier mesons. An ?7/7r -ratio of 0.55 ± 0.08 was as- 
sumed @ , consistent with the measured result presented 
here. 
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FIG. 3: The ratio of the r\ and n° yield as a function of pt 
in p+C collisions at ^/snn = 17.4 GeV. The 7r° spectrum 
is taken from [TJ|. The error bars represent the quadratic 
sum of the statistical and systematic uncertainties. The ratio 
assumed in [t| with its respective uncertainties is indicated 
by the horizontal lines. 



The ^-production in p+Pb and p+C collisions at 
y/snn = 17-4 GeV has been analyzed via the 77 — > 2-f- 
channel by an invariant mass analysis analogous to that 
used for the neutral pions [l8T | : the two-photon invariant 
mass is calculated for all possible photon pair combina- 
tions in an event, and the r\ yield is extracted statistically 
from the counts in the 77 mass peak. A requirement that 
the asymmetry of the photon pairs a = , with E\^, 

being the energies of photon 1 and 2, should be less than 
0.7 was applied to reduce the combinatorial background. 
In order to remove this background, the invariant mass 
distributions of uncorrclatcd photon pairs taken from dif- 
ferent events with identical multiplicity were determined. 
These distributions were scaled to the background level 
in the real pair distribution for each px-mterval and sub- 
tracted to obtain the pr-dependent ri-yield. The accep- 
tance and efficiency corrections to the rj yield were deter- 
mined from the same Monte- Carlo simulations used for 
the inclusive photon and neutral pion analyses. 

A significant r\ result could be obtained only for the 
p+C data set with the HEP trigger. The fully corrected 
r\ spectrum is shown in Fig. [2] The resulting measured 
?7/7r -ratio for p+C shown in Fig. [3] is consistent with the 
ratio assumed in [Sjl, shown as horizontal lines. 
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FIG. 4: Double ratio i? 7 = {j/-K ) me as/{j/-K°)decay for p+Pb 
(upper panel) and p+C (lower panel) collisions at ^/snn = 
17.4 GeV. A fit to the neutral pion data from [JJ| has been 
used for the measured tv° result. The error bars show the 
statistical errors and the shaded areas show the systematic 
uncertainties of the measurements. 



The double ratios i? 7 for p+Pb and p+C are shown 
in Fig. |U The uncertainty of the energy scale largely 
cancels out in the double ratio. No significant photon 
excess is observed in the pt range relevant for thermal 
photon production. However, it is possible to determine 
upper limits for the production of direct photons up to 
Pt = 3.2 GeV/c, which are shown in Fig. [SJ The upper 
limits are extracted as 1.28c above the measured result 
(using i? 7 = max[l,i?™ eas lh as used for the previously 
published Pb+Pb results [|, i, and are plotted as 
vertical arrows with the tail indicating the upper limit. 

These new results are compared to the previously pub- 
lished results for central Pb+Pb collisions @, Q in Fig. [6] 
after scaling the upper limits for the p+A datasets by 
the corresponding ratio of the number of binary nucleon- 
nucleon collisions, (N co u), as determined from Glauber 
calculations [20| using a nucleon-nucleon inelastic inter- 
action cross section of er^^ = 31.8 mb. The transverse 
energy was modeled by sampling a negative binomial 
distribution to determine the Et contribution of each 
participating nucleon. To these the same selection on 
Et as used in the event selection by the WA98 trigger 
was applied, to determine the (N co u) values. With this 
procedure the number of binary nucleon-nucleon colli- 



sions has been determined as (N c 
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FIG. 5: Upper limits, quoted at 1.28cr, for the direct photon 
production in p+Pb and p+C collisions at ^/%n = 17.4 GeV. 
The p+C values have been shifted higher by 0.05 GeV/c in 
Pt for visibility. 
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TABLE I: Systematic uncertainties on the direct photon yield 
measurement at two representative px values. 



(N coU ) p+Pb = 3.8 + 0.4 [18]. For the Pb+Pb data, the 
events were selected as the 13% most central fraction of 
the minimum bias cross section of o ,p ^ rPb ~ 6300 mb, 
corresponding to (N coU ) Pb+Pb ^ cent = 727.8 + 72.8. 

The direct photon yields observed for central Pb+Pb 
collisions, scaled by the number of nucleon-nucleon colli- 
sions, are consistent with the upper limits from both the 
p+Pb and p+C results, as shown in Fig. [6l Therefore, 
the new results do not allow to constrain the prompt pho- 
ton contribution to the Pb+Pb direct photon yield, and 
hence do not allow to draw definitive conclusions about a 
possible thermal contribution to the Pb+Pb direct pho- 
ton yield. A new NLO pQCD calculation is shown in 
both panels for p+p collisions at ^/snn = 17.4 GeV [2lj . 
scaled by (N co u) Pb+Pb cent . The calculation employs the 
GRV photon fragmentation function [22| and is shown 
for scales fi — 0.5pr,Pr, %Pt, which reflect the theoreti- 
cal uncertainties. Within these the NLO calculations are 
able to describe the central Pb+Pb direct photon excess, 
and hence do not provide evidence for a thermal direct 
photon contribution. 




p T (GeV/c) 



FIG. 6: Comparison of the direct photon spectra for central 
Pb+Pb collision from [9j to the p+Pb (upper panel) and p+C 
results (lower panel) presented in this paper. The two p+A 
data sets have been scaled by the ratio of the number of binary 
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lations for p+p collisions scaled by N^ b f b for three different 
scales are shown as lines. 



In summary, new measurements of r\ production in 
p+C collisions at ^/snn = 17.4 GeV have been presented. 
The result is in agreement with previous ?7/7r -ratio re- 
sults. Inclusive photons have been measured in p+Pb 
and p+C collisions at ^/snn = 17.4 GeV, and upper lim- 
its on the production of direct photons were obtained. 
The results cover the pr range of 0.7 GeV < pt < 
3.2 GeV/c. Comparison of these upper limits, scaled 
by (N co u), to the results obtained by WA98 for central 
Pb+Pb collisions at v /snn = 17.3 GeV does not improve 
the existing experimental constraints on the production 
of thermal direct photons at SPS energies. Similarly, 
new NLO pQCD calculations of the prompt photon con- 
tribution are consistent, within uncertainties, with the 
observed central Pb+Pb direct photon excess. Hence, it 
remains an open question to what extent the direct pho- 
ton signal observed by WA98 for central Pb+Pb colli- 
sions at i/snn = 17.3 GeV may be attributed to thermal 
radiation. 
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